"t^" ■ Long- Range Correlation of Electron Pairs in the Hubbard 



0^ 






o 



Model at Finite Temperatures in Three Dimensions 



O, 
^ '_ Gang Su and Bao-Heng Zhao 

Department of Physics, Graduate School 

>• ■ Chinese Academy of Sciences, P. O. Box 3908 

2 ; Beijing 100039, China 



Abstract 



We show that in the translation invariant case and in the antifer- 



'^ ' romagnetic phase, the reduced density matrix p2 has no off-diagonal 

O ■ long-range order of on-site electron pairs for the single-band Hubbard 



model on a cubic lattice away from half filling at finite temperatures 



^ . both for the positive coupling and for the negative coupling. In these 

H ' 

CZ \ cases the model can not give a mechanism for the superconductivity 

caused by the condensation of on-site electron pairs and the nearest- 
neighbor electron pairs. 

PACS numbers: 71.10.-Fx, 74.20.-z, 05.30.-d. 



It was proposed the possibility that the Hubbard model can give a mech- 
anism to characterize high temperature superconductivity [Q. It has been 
shown that the rj pairing states for the single band Hubbard model with con- 
stant hopping matrix element for nearest neighbor sites have ODLRO|0]. It 
has also been shown that in the ground state the model with negative cou- 
pling has off-diagonal long-range order (ODLRO) on some bipartite lattice^]. 
However, as the temperature T > whether the reduced density matrix p2 
has ODLRO is inconclusive. The purpose of this work is to study the ODLRO 
of p2 in the single-band Hubbard model with constant hopping matrix ele- 
ment for nearest neighbor pairs on a cubic lattice at finite temperatures. In 
the translation invariant case we obtain an equation ( (10) in the following) 
which is satisfied by the long-range correlation functions of on-site electron 
pairs. From which we find that there is no ODLRO of on-site electron pairs. 
Then we consider the case with an antiferromagnetic background. In this 
case we will show that there is still no ODLRO of on-site electron pairs. 

The Hamiltonian of the model is 

H = -t Y. {ala^> + h%) + UY.ala,h% - /i^(a|;ar + h%), (1) 

(r,r') 

where Or (&r) is the annihilation operator of the electron with spin up (down) 
at site r. U and t are constant, and /i is the chemical potential, (r, r \ is a 
pair of nearest-neighbor sites. The lattice is cubic in the three dimensional 
space. We take the lattice spacing to be unity. 
Define the ?]— operators as those in [@, §, ^, 

^- = E e"'""ar&r, V+ = ^-, Vz = - E(4«r + h% - 1), (2) 



where 7r= (vr, 7r,7r). They satisfy 

[V+,V-] = '^Vz, [v±,Vz] = TV±, (3) 

and 

r]^H ={H + U-2i2)r]^. (4) 

Consider the average over a grand canonical ensemble 

iv+V-) = ^Tr(77+r7_e-^^), (5) 

where Z is the grand partition function; [3 the inverse temperature. Using 
(4) and cyclically permuting the factors under the trace, we obtain {Ti-r]^) = 
{r]^r]_) e^^^~'^^\ Then the first one of (3) gives, when U ^ 2/i, 



From (2) we have (r^^) = {N — M)/2, where N = J2r (ajflr + ^J^r) is the aver- 
age total number of electrons, M the total number of lattice sites. Denoting 
the density of the number of electrons as -D = N/M, (6) leads to 

f/ 7^ 2/i corresponds to the case away from half filling. We will only consider 
this case below. As T = /3~^ = finite and away from half filling, the right 
hand side of (7) shows that -p {Tj^rjJ) is a finite intensive quantity, which is 
still finite in the thermodynamic limit ( M ^ oo, A^ ^ oo with D = N/M 
and f3 fixed). 

On the other hand, by the definition of ri±, we can write 



First we consider the case without the spontaneous breaking of the transla- 
tion symmetry on the lattice, and take the periodic boundary condition. We 
have 

j^{v+V^) = T.{bWr<^obo)e'^-. (9) 

In the thermodynamic limit, the right hand side of (9) becomes an infinite 
series. As we show above that -^ (v+V-) is finite, thus the infinite series 
J2\blalaobo) e*'^''" is convergent . 

Let us take the thermodynamic limit first, then, for the aim to investigate 
the behavior of (blalaobo) at large |r| , consider a cube of volume (25*)^ in 
the infinitely large system. The center of the cube is at the origin of the 
coordinates. Let f{L) = J2 (blalaobo) e^^''^ , where r = {l,m,n), —L < 

faces 

l,m,n < L, 1,171,71, L are integers, and L < S. X) is fhe sum over the six 

faces 

faces of the cube with volume (2L)^, i.e., 

L L L 

E= E + E + E • 

faces l,m=—L{n=±L) n,l=—L(m=±L) m,n=—L{l=±L) 

It is clear that the sum of (blalaobo) e'^''^ over the cube of volume (25')'^ can 
be written as J2 (blalaobo) e'^''^ = J2 f{L)- Since as S* ^ oo this series is 

r \ ' L=0 

convergent, we have limi„^oo/(-^) = 0, namely, 

lim V (^iajaofoo) e*"-' = 0. (10) 

faces 

Eq.(lO) is a strong restriction on the long-range correlation functions. We 
will show that it determines (blalaobo) completely in the limit |r| -^ oo. 

From the algebraic approach to the quantum statistical mechanics we 
know that the correlation functions in the equilibrium states of a pure ther- 
modynamic phase have the spatial cluster properties , which is a rigorous 



result IHI. It enables us to write 

[blalaobo) -^ (&J4) {(^oK) , as \r\ -^ oo. (11) 



In the translation invariant case (flr^r) = {o-obo)-, thus ih\a\aoho) -^ \{0'obo)\ , 
as |r| -^ oo. In (11) (or^r) should be understood as Bogolubov's quasi- 
average, i.e., a f/(l) symmetry breaking term was added to the Hamilto- 
nian, and after taking the thermodynamic limit it has been sent to zero0]. 
{aj.br) does not vanish, if there is on-site electron pair condensition, other- 
wise it vanishes. In the calculation of (blalaobo) the symmetry breaking term 
is not necessary. Since blalaobo is an invariant under the global U{1) gauge 
transformation, its usual average ( in the Hamiltonian without the symmetry 
breaking term) and the quasi- average are the same0. 

It is a well established empirical fact that in the second-order phase tran- 
sition the correlation length is divergent at the critical point, and it is finite 
away from the critical point. The modern theories of critical phenomena (the 
scaling hypothesis, the renormalization group approach, etc.) are based on 
it|^. We assume that this empirical fact works in the present problem, which 
is consistent with the superconducting transition being a second-order phase 
transition. Define G(r) = (blalaobo) — \{aobo)\ for arbitrary r. Suppose that 
the system is in a noncritical thermodynamic state, then G{r) can not decay 
slower than 0(|r|~ ) as |r| — > oo, otherwise the correlation length for G{r) 
will diverge to infinity at a noncritical point. So limi_>oo J2 G{r)e'^'^''' = 0. 

faces 

Hence (11) shows that we can write 

lim Y. (blalaobo) e'''-''=\{aobo)f lim Y. e'""- (12) 

L— >oo r — \ ' L^oo ' — 

faces faces 



It can be shown that 

E e''^" = (-1)^2, for L > 1. (13) 

faces 

Then (10) leads to (—1)^2 |(ao&o)| -^ 0. Therefore we obtain {aobo) = 0, and 
by translation invariance, 

{ttrbr) = 0, at any site r. (14) 

Furthermore (11) and (14) show that 

blalaobo) — > 0, as |r| — > oo. (15) 



Eq.(14) shows that there is no condensation of on-site electron pairs, and 
(15) shows that the reduced density matrix p2 has no ODLRO for on-site 
electron pairs. These two statements are equivalent. 

Using (14) and taking the quasi-average of [H, aj.br] , it can be shown 

(arV) = 0, (16) 

where r, r are nearest neighbor sites. From the cluster property of correlation 



or 1, 



s — s 



functions, we can write (bla^agb^') — > 0, where 
finite, and |r — s| ^ oo. 

The above results are obtained for the translation invariant case. The an- 
tiferromagnetic order breaks the translation symmetry. But we can show that 
the existence of antiferromagnetic order, if any, is not in contradiction with 
(14)-(16). In the antiferromagnetic phase there are two kinds of translations, 
i.e., the displacement Ar = (Ax, Ay, Az), Ax + Ay + Az = even integer. 
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and Ax + Ay + Az = odd integer. The first one does not break the transla- 
tion symmetry, but the second one does. The second kind of translations is 
equivalent to exchange the up spins and the down spins, so it is easily seen 
that 

(ar6r) = e*""" (aobo) , (blalasK) = (bl-salsaobo) ■ (17) 

Thus (9) and (10) are still valid, but (12) should replaced by 

limi^oo E (blalaoboje^''''' = \{aobo)f\imL^oo E 1 

faces faces fl8l 

= |(aA)|'limL^oo2(12L2 + l). 

So we can show that (14)-(16) still hold, i.e., there is no ODLRO for on-site 
electron pairs. 

In the derivation of (14)-(16), the sign of U plays no role whatsoever, so 
the above results apply to the cases both for f/ > and for U < 0. 

In summary, we have shown that for the translation invariant case or in 
the antiferromagnetic phase the reduced density matrix p2 has no ODLRO 
of on-site electron pairs for the single band Hubbard model on a cubic lattice 
away from half filling at finite temperatures both for f/ > and for f/ < 0. 
In these cases the model can not give a mechanism for the superconductivity 
caused by the condensation of on-site electron pairs and nearest neighbor 
electron pairs. To derive the above results we have assumed that in the case 
away from the critical point the correlation length is finite , and no other ad 
hoc assumption is needed. The above discussion can be easily generalized to 
any dimensions. Our results are not incompatible with f^ and |^ , since their 
results are not the ensemble average at T > 0. 
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